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Abstract —The main challenge of ultra-reliahle machine-to- 
machine (M2M) control applications is to meet the stringent 
timing and rellahllity requirements of control systems, despite 
the adverse properties of wireless communication for delay 
and packet errors, and limited battery resources of the sensor 
nodes. Since the transmission delay and energy consumption of a 
sensor node are determined by the transmission power and rate 
of that sensor node and the concurrently transmitting nodes, 
the transmission schedule should be optimized jointly with the 
transmission power and rate of the sensor nodes. Previously, 
it has been shown that the optimization of power control 
and rate adaptation for each node subset can be separately 
formulated, solved and then used in the scheduling algorithm in 
the optimal solution of the joint optimization of power control, 
rate adaptation and scheduling problem. However, the power 
control and rate adaptation problem has been only formulated 
and solved for continuous rate transmission model, in which 
Shannon’s capacity formulation for an Additive White Gaussian 
Noise (AWGN) wireless channel is used in the calculation of the 
maximum achievable rate as a function of Signal-to-Interference- 
plus-Noise Ratio (SINK). In this paper, we formulate the power 
control and rate adaptation problem with the objective of 
minimizing the time required for the concurrent transmission 
of a set of sensor nodes while satisfying their transmission 
delay, reliability and energy consumption requirements based 
on the more realistic discrete rate transmission model, in which 
only a finite set of transmit rates are supported. We propose a 
polynomial time algorithm to solve this problem and prove the 
optimality of the proposed algorithm. We then combine it with 
the previously proposed scheduling algorithms and demonstrate 
its close to optimal performance via extensive simulations. 

Index Terms —ultra-reliable communication, machine-to- 
machine communication, power control, rate adaptation, energy 
efficiency, delay constraint. 

I. Introduction 

M2M is a new communication paradigm that aims to enable 
the communication between the devices with no or mini¬ 
mal human interaction. The wireless communication between 
sensor nodes and controllers in ultra-reliable M2M control 
applications reduces the cost of their installation, maintenance 
and upgrade, together with the complexity of the overall 
system, compared to their wired equivalent m. The main 
challenge of ultra-reliable M2M control applications is the 
design of a robust scheduling algorithm that satisfies the 
packet generation period, transmission delay and reliability re¬ 
quirement of the sensors required to maintain a certain control 
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system performance 0, B despite the adverse properties of 
wireless communication for non-zero packet error probability 
and non-zero delay at all times. Since the transmission delay 
and energy consumption of a sensor node are determined by 
the transmission power and rate of that sensor node and by the 
same of the other nodes that are scheduled for concurrently 
transmission, the schedule should be optimized jointly with 
the transmission power and rate of the sensor nodes. 

Ultra-reliable communication (URC) is defined as a high- 
level communication service which is available almost 100% 
of the time. a, Q discuss the relation between control 
information and actual data in URC based systems. pro¬ 
poses a URC based system based on the prediction of future 
reliable channels. None of these works however consider the 
periodic packet transmission of the sensor nodes in control 
applications. 

Joint optimization of power control, rate adaptation and 
scheduling has been studied for both general purpose wireless 
networks Q, ii, 0, cni and wireless networked control 
systems (WNCSs) HI], ifT^ . This optimization problem has 
been formulated and solved for general purpose networks with 
the goal of satisfying either a single deadline for all packets 
Q, in or individual deadlines for each packet 0, lITo] . In 
addition to the delay constraints, the formulation of this joint 
optimization problem for WNCSs has included the uniform 
distribution of packet transmissions over time to satisfy the ro¬ 
bustness requirement of the control systems, and the periodic 
data generation, reliability and energy consumption require¬ 
ments of the sensor nodes ifTTl . ifT^ . It has been shown in ifT^ 
that the optimization of power control and rate adaptation for 
each node subset can be separately formulated, solved and 
then used in the scheduling algorithm in the optimal solution 
of the joint optimization of power control, rate adaptation 
and scheduling problem. However, the power control and rate 
adaptation problem has been formulated for continuous rate 
transmission model uni, HD, in which Shannon’s channel 
capacity formulation for an AWGN wireless channel is used in 
the calculation of the maximum achievable rate as a function 
of SINR. In practical communication systems, however, only 
a finite set of discrete transmission rate values are supported, 
as suggested by the practical realization of multiple data rates 
in ifTSll . 

The goal of this paper is to study the power control and 
rate adaptation problem with the objective of minimizing 
the time required for the concurrent transmission of a set 
of sensor nodes while satisfying their transmission delay, 
reliability and energy consumption requirements based on the 
more realistic discrete rate transmission model. We propose 
a polynomial time algorithm for the power control and rate 


adaptation problem and provide the proof of the optimality of 
the proposed algorithm. We then demonstrate the performance 
of the proposed algorithm within the previously proposed 
scheduling framework for ultra-reliable M2M control via 
extensive simulations. 

The rest of this paper is organized as follows. Section HI] 
describes the system model and assumptions used throughout 
this paper. Section uni formulates the power control and rate 
adaptation problem, proposes a polynomial time algorithm and 
proves the optimality of the proposed algorithm. Performance 
results of the proposed algorithm are presented in Section HVl 
The hnal review of our work and the ideas for future research 
are given in Section fVl 

11. System Model and Assumptions 

1) The network contains sensor nodes, controllers and 
actuators. The controllers can receive packets only from 
one sensor node at a time. When a controller receives 
the sensor data, it performs the new control computa¬ 
tion and sends the output to the actuator. We assume 
that the actuators always receive the controller packets 
successfully. 

2) We assume that sensor nodes transmit their data pack¬ 
ets to their corresponding controllers in one hop. The 
extension for the multi-hop packet transmission of the 
packets from the sensor nodes to the controllers is out 
of scope of this paper and subject to future work. 

3) Central controller is selected from the existing con¬ 
trollers and responsible for network synchronization, 
resource allocation and scheduling of the sensor nodes. 
Central controller is assumed to have complete network 
topology, however the method for topology learning is 
out of scope of this paper. 

4) TDMA is used as a Medium Access Control (MAC) 
protocol due to its superior delay and energy perfor¬ 
mance for the networks with predetermined topology 
and data generation pattern m. The time is divided 
into subframes, which are further divided into beacon 
and time slots. The central controller sends the beacon 
at the beginning of each subframe and includes the 
scheduling and resource allocation information of each 
scheduled node within the subframe. If there exist 
sensor nodes outside the transmission range of the 
central controller then the central controller uses the 
corresponding controllers to relay the beacon to these 
nodes. 

5) We only consider the energy consumption for the data 
packet transmission of each sensor node since the 
energy consumption in active mode is much larger 
than that in sleep and transient modes, and the energy 
consumption in the reception of the beacon packets is 
not subject to optimization. We define the maximum 
allowed per packet energy consumption to achieve a 
certain lifetime by e/ for sensor node 1. 

6 ) The transmission power can take any value below Pmax- 
This continuous power assumption is frequently used in 
previous scheduling algorithm designs since practical 
radios support a large number of discrete power levels. 


7) 


8 ) 


9) 


resulting in high approximation accuracy with lower 
complexity. 

We assume that the channel is perfectly known and 
constant over the scheduling frame, i.e. channel gain gik 
from the sensor node I to the controller corresponding to 
sensor node k does not change. Considering the channel 
variation of the links within the scheduling frame so fast 
fading wireless channels is out of scope of this paper 
and subject to future work. 

We use the discrete rate transmission model in which a 


hnite number of transmission rates r = (r^,r^ 


rQ) 

corresponding to a hnite number of SINR levels 7 = 
are determined such that node I can 
transmit at rate if the SINR achieved at the corre¬ 
sponding link, i.e. 7 / = No+j:i^lpigik ’ greater than 
or equal to 7 ^, where pi is the transmit power of node I, 
Nq is the background noise power and Q is the number 
of discrete transmission rate levels. Each rate level r®, 
i € [1,Q] is determined based on SINR-rate mapping 
function /( 7 ), i.e. r® = /( 7 *), i S [1, Q], which satishes 
the constraint given by ^ < 0. This constraint is 

satished for Shannon’s capacity formulation which is 
commonly used in AWGN channels. 

The scheduling algorithms proposed for WNCS in IIT2II 
has been adopted to be jointly used with the proposed 
power control and rate adaptation algorithm. The goal 
of the scheduling algorithm is to provide maximum 
robustness accommodating packet losses, topology and 
sensor requirement changes. The subframe and frame 
length are defined as the minimum and maximum packet 
generation period of all sensor nodes, and denoted by 
S and F, respectively. The total active length of a 
subframe m, denoted by Om, is defined as the sum 
of the length of the time slots allocated to subframe 
m. The robustness of the schedule requires distribution 
of data transmissions as uniformly as possible over the 
scheduling frame. The objective is therefore quantified 
as minimizing the maximum total active length of all 
the subframes in a frame. 

To illustrate the scheduling, let us assume the example 
network consist of 4 sensor nodes. Sensor node 1 has 
packet generation period of 1 ms and packet transmis¬ 
sion time of ti = 0.15 ms. Sensor nodes 2, 3 and 4 
have packet generation period of 2 ms and their packet 
transmissions require time slot allocations of lengths 
t 2 = 0.20 ms, fa = 0.25 ms and = 0.30 ms, 
respectively, when no concurrent transmissions are al¬ 
lowed. Let us assume that the concurrent transmission of 
nodes 2 and 3 takes 0.30 ms, i.e. ^ 2,3 = 0.30 ms. Here, 
the subframe length 5” = 1 ms and the frame length 
F = 2 ms, containing 2 subframes. Fig. [T|shows the ro¬ 
bust schedule that minimizes the maximum total active 
length of the subframes in the presence of the concurrent 
transmissions. The maximum total active length is equal 
to max(ai,a 2 ) = inax{ti -f ^ 2 . 3 , ti + Ia) = 0.45 ms. 




-ai-► - 32 - 

-ti--12,3-► -14- 


Sensors 2,3 Sensor 1 Sensor 4 


-Subframe 1 (1 ms)--Subframe 2 (1 ms)-► 

Fig. 1. Illustration of the robust scheduling algorithm for the network 
of 4 nodes. 

The scheduling algorithms proposed in IfT^ are based 
on first assigning the nodes to the subframes through 
the node assignment algorithm and then determining 
the subsets of nodes that concurrently transmit through 
the concurrency allocation algorithm. The first algorithm 
adopts Sorted Node Assignment (SNA) and Minimum 
Length Allocation (MLA) hence is denoted by SNA- 
MLA algorithm, whereas the second algorithm adopts 
SNA and Maximum Utility Allocation (MUA) therefore 
is called SNA-MUA. SNA is based on assigning each 
node to the subframe with the minimum total active 
length giving priority to the nodes with higher trans¬ 
mission times and considering the nodes connected to 
different controllers separately. MLA algorithm enumer¬ 
ates all feasible subsets of the nodes with the same 
packet generation period assigned to the same subframe, 
and then selects some of these feasible subsets such that 
every node is covered at least once and the total time 
allocated is minimized. MUA algorithm, on the other 
hand, iteratively constructs the concurrently transmitting 
node subsets from the nodes with the same packet 
generation period assigned to the same subframe by 
including the node that maximizes the utility, which is 
defined as the decrease in the transmission time of a 
set of sensor nodes when they transmit concurrently, in 
the subset. The SNA algorithm uses power control and 
rate adaptation in determining the transmission time of 
each sensor node while satisfying its delay and energy 
requirement whereas the MLA and MUA algorithms use 
power control and rate adaptation algorithm in deciding 
whether the concurrent transmission of a subset of sen¬ 
sor nodes while satisfying their transmission delay and 
energy requirements is feasible and their transmission 
time if feasible. 

III. Power Control and Rate Adaptation Problem 
A. Problem Formulation 

The power control and rate adaptation problem aims to 
minimize the time required for the concurrent transmission 
of the sensor nodes, while satisfying their transmission delay 
and energy consumption requirements. As described in detail 
in Section m this problem needs to be solved to provide 
the transmission time of a single node and subset of nodes 
in the node assignment and concurrency allocation parts, 
respectively, of the SNA-MLA and SNA-MUA algorithms. 
The formulation of the power control and rate adaptation 
problem for the sensor nodes in the set S is given as 
minimize 

t (la) 
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(Ic) 

tipi < ei, 

1 £ S 

(Id) 

pi S Pmax , 

1 £ S 

(le) 

, Ri 
ti = —, 

Xl 

1 £ S 

(If) 

kiqy'^{No + '^Pugui) > 0, l£S,q£[l,Q] 

U^l 

(Ig) 

Q 

9=1 

1 £ S 

(Ih) 

Q 

^ klqp 

l£S 

(li) 


9=1 


variables 

t > 0, Pi > 0, kiq e {0, 1}, i e s, q £ [i, Q] 

where t is the time duration required for the concurrent 
transmission of the nodes in the set S', ti and xi represent 
the transmission time and the transmission rate of node I, 
respectively; Ri is the length of the packet of sensor node (; 
kiq is the indicator of the assignment of sensor node I to the 
q-th rate level; di is the delay requirement of sensor node I 
to achieve a certain control performance. Eqs. ([Tal l and (Sbli 
together represent the objective of minimizing the time slot 
length allocated for the concurrent transmission of the nodes, 
i.e. minimizing maximum transmission time of the nodes. Eqs. 
(fT^ and (fTdl i represent the delay and energy constraints of the 
nodes, respectively. Eq. (fTel i states that the transmit power of 
the Lth node cannot be greater than the maximum allowed 
power. Einally, Eqs. (fTft - (|Ti|l represent the computation of the 
transmission time of node I by assigning one of the discrete 
transmission rates in the set r that satisfy the corresponding 
SINR requirement. 

Next, we propose a polynomial time algorithm for the 
solution of this optimization problem and give a detailed proof 
of the optimality of the proposed algorithm. 

B. Proposed Polynomial Algorithm 

We propose Longest Transmission Time Eirst (LTTE) al¬ 
gorithm to determine the optimal power and rate allocation 
for the concurrent transmission of the nodes in the set S in 
polynomial time. The algorithm is based on the initialization 
with the vector of the lowest possible rates satisfying the 
delay constraints of the nodes and then checking the feasibility 
of the lowest possible number of rate allocation vectors for 
satisfying energy, maximum power and SINR constraints. 

Let us enumerate the links in the set S from 1 to [S'!. Let 
r be a vector of the transmission rates of the links in the 
set S. Determining the feasibility of a rate vector r requires 
checking the existence of a power vector corresponding to r, 
delay, energy and maximum power constraints. Since the rates 
of the links are predetermined by the corresponding elements 


ti<t, I £ S 















of the vector r, the corresponding SINR threshold levels are 
predetermined. Then the existence of a transmit power vector 
that satisfies the constraint 

Pi9ii 

No + Y.j^^Po9j^ 

for every node i, where 7 ^ is the SINR level corresponding to 
ri, can be found by testing Perron-Frobenius conditions M- 
Perron-Frobenius conditions determine whether there exist a 
power vector satisfying these SINR constraints, and if so, the 
component-wise minimum power vector p"*®**. is then 

tested for delay, energy and maximum power constraints given 
by Eqs. (fT^ . (fTdl i and ([Tel l, respectively. 

Let t denote the time duration required for the concurrent 
transmission of the nodes in the set S. The optimal values of 
r and t are denoted by and respectively. The i-th 
element of the vector corresponding to the optimal rate 
of node z, denoted by r°^*, is initialized to the minimum rate 
level that satisfies the delay constraint of node i, denoted by 
given by Eq. (fT^ for all z G [ 1 , 1 'S'I], whereas is 
initialized to infinity (Lines 1 — 2). If this initial is not 
feasible, then the algorithm returns = 00 , meaning that 
there does not exist a feasible solution (Lines 3 and 15). Oth¬ 
erwise, the algorithm calculates the time duration of the time 
slot corresponding to the rate vector by determining the 
maximum transmission time of the nodes in the set S (Lines 
4 — 5). Since the node with the maximum transmission time 
denoted by j determines the time slot length, the length of 
the time slot can only decrease by increasing its transmission 
rate. Therefore, the algorithm continues by increasing the rate 
of the j-th node by one level if it is less than the maximum 
possible rate (Lines 6—7), and then checking the feasibility 
of the resulting vector (Line 3). The algorithm stops if 
either the j-th node has rate meaning that it is not possible 
to decrease the time slot length any further (Lines 8 — 9), or 
the resulting vector is not feasible. 

The complexity of the LTTE algorithm is 0{Q\S\^) since 
the algorithm checks at most Q|S'| transmission rate vectors 
for feasibility and the complexity of the feasibility check 
is dominated by the complexity of testing Perron-Lrobenius 
conditions of complexity 0(15^) ifTSll . 

C. Optimality of Proposed Algorithm 

We now prove the optimality of the proposed algorithm 
following the statement of the following three lemmas. 

Let r = (ri, r 2 ,..., r|s|) be a vector of the transmission 
rates of the links in the set S and t be the length of the 
time slot corresponding to the rate vector r, such that t = 
TOaa;ig[i | 5 |]-p^, if r is feasible. Denote the rate index vector 
j.d _ such that for at least one element j, 

Tj < Tj, while rf < ri for all i G [ 1 , IS”!], by the descendant 
of r. 

Lemma 1. Let both r and r'^ satisfy the delay constraint given 
in Eq. (O- If r'^ is infeasible, then r is infeasible. 

Proof: Let us denote the component-wise minimum 
power vectors corresponding to the rate vectors r and by 
P = (Pi,P 2 ,---,P|S|) and p^ = {pi,pi,...,pfg^), respectively. 


Algorithm 1 Longest Transmission Time Lirst (LTTE) algo¬ 
rithm_ 

Input: Node set S considered for concurrent transmission 
Output: Optimal rate values for the nodes in the set S if their 
concurrent transmission is feasible 


1 : 

2 : 

3: 

4: 

5: 

6 : 

7: 

8 : 

9: 

10 : 

11 : 

12 : 

13: 

14: 

15: 


= 00 ; 

_ {^th ^th ^th \ 

r _ , ^2 , 

where rf = TOm,g[i,Q]{r9|^ < dj; 
while is feasible do 

j = arp7TzaXig[i,|s|]^; 



increase by one level; 
else 
break; 

end if 
end while 

if < 00 and r°P* is not feasible then 
decrease rf* by one level; 

end if 
return r°P^, 


We first show that pf < pi for all z G [1, l^j]. Denote the 
SINR threshold levels corresponding to the rates and rf 
by and yf, respectively. Since higher rate corresponds 
to a higher SINR threshold value and rf < ri, yf < 7 ^. 
Therefore, any power vector satisfying the SINR requirements 
for the SINR threshold levels ^,1 G [1, liSI] also satisfies 
these requirements for yi,i G [1, l^l]. Suppose that pf > 
for an arbitrary link k, then p'^ cannot be component-wise 
minimum power vector. This is a contradiction. 

If is infeasible due to the maximum power constraint 
given in Eq. (HHl, then there exists at least one link k for 
which pf > Pmax- Since pk > pf, Pk > Pmax- Hence, r does 
not satisfy maximum power constraint either. 

If r'^ is infeasible due to the energy constraint given in Eq. 
dMll, then there exists at least one link k for which tfpf > e^. 
Let a denote the ratio of and rf. Since SINR-rate mapping 
function is an increasing function with a non-negative second 
derivative as described in Section HIl the transmission power 
Pk will be greater than apf. Therefore, tkPk > -^otpf = 
tfpf > Cfc so r is infeasible due to the energy constraint. 

■ 

Lemma 2. Let and t are the transmission times corre¬ 
sponding to the rate index vectors and r, respectively. Then, 
>t is always true. 

Proof: We know that, t = rnaxi^[i^is\]ti = 
Since we are considering the case where 
rf < ri for all z G [1, IS”!], we can simply show that t = 
TOaa;ig[i,|s|]-^ < waa:ig[ijs|]^ Hence, P^ >t. ■ 

Lemma 3. Suppose that node j has the maximum time slot 
length among concurrently transmitting [S'! nodes, i.e. j = 
argmaXi^[i^\Swf^- Then, increasing the transmission rate of 
any node k f j, where k,j G [1, [S'!], does not decrease the 








time slot length t. 

Proof: t = ^ ^ \/k f j, where 

k,j G [1,1 S'!]. Suppose that, for an arbitrary link k 7 ^ j, we 
increase to Then -% < -^ < = t. Hence, 

Tk Tk — rj 

increasing the transmission rate of any node k j does not 
decrease t. ■ 

Theorem 1. Algorithm stated in Section IV gives the optimal 
result. 

Proof: Suppose that LTTF returns the as an optimal 
solution and the j-th link has the maximum transmit time 
among all the sensor nodes. Now, consider that is not an 
optimal solution. This means that there exists another optimal 
rate index vector r** such that t** < t°P*. 

Due to Lemma 2, the transmission times corresponding 
to all the descendants of r°^’‘ are greater or equal to that 
corresponding to r°P*^ itself, i.e. t‘^ > t°P*. This means that 
r** cannot be a descendant of r°P*. Also due to Lemma 3, 
increasing the rate value of any other node than j cannot 
decrease 

We now introduce the new rate index vector in 

which only rj is increased by one level, i.e. = 

r 2 ^*,..., ..., Since LTTF returns we 

know that r°P‘’+ is infeasible. Due to Lemma 1, for at least 
one sensor k 7 ^ j, r^* < must hold for r** to be feasible. 
We also know that at each iteration of the LTTF algorithm, the 
rate value of the sensor node with the maximum transmission 
time is increased by one level. Therefore at some iteration of 
the algorithm, the fc-th node is evaluated as k = argmax^ 
Vi € [1,15*1] and rf* increased by one level. It follows that, 
t** > = t°P*, which is a contradiction. 

■ 

IV. Simulation Results 

Simulations are performed in MATLAB. The sensor nodes 
and the controller are uniformly distributed within a square 
area. The figures show the average of 100 independent random 
network topologies. The packet generation period and the 
packet transmission requirement of each node are randomly 
chosen from the sets [1,2,4, 8 ] and [50,100], respectively. 
Number of controllers is fixed and equal to 3. The continuous 
rate model adopts Shannon’s channel capacity formulation for 
an AWGN wireless channel in the calculation of the maximum 
achievable rate as a function of SINR and is denoted by cont 
in the legends of the figures. Discrete rate model adopts 4 
discrete rate levels corresponding to 4 SINR levels equal to 
[—c», 10,20,30] dB and 8 discrete rate levels corresponding 
to 8 SINR levels equal to [—c», 0, 5,10,15, 20,25, 30] dB cal¬ 
culated by using Shannon’s capacity formulation, and denoted 
by disc4 and discS, respectively, in the legends of the figures. 
The channel attenuation is calculated by using Rayleigh fading 
with scale parameter set to the mean power level determined 
by using the large scale statistics modeled as 

PL{d) = PL{do) + 10a \og{d/do) + Z, (3) 

where d is the distance between the node and the controller, 
PL{d) is the path loss at distance d, PL{df) is the path loss 


at reference distance c?o, Z is the Gaussian random variable 
with zero mean and standard deviation cr^. All the simulation 
parameters are given in Table U 

TABLE I 

Simulation Parameters 


PL{do) 

70 dB 

Pmax 

250 mW 

a 

3.5 

O-z 

4 dB 

No 

10"® W/Hz 

w 

100 MHz 


Fig. |2] shows the normalized maximum total active length of 
SNA-MLA and SNA-MUA scheduling algorithms for continu¬ 
ous and discrete rate transmission models for different number 
of nodes at 5 nodes/m^ density. The normalized maximum 
total active length is defined as the ratio of the maximum 
total active length of the scheduling algorithm to that of the 
optimal solution for the continuous rate transmission model. 
The reason for normalizing the maximum total active length 
of even the discrete rate by that of the continuous rate 
transmission is to observe the performance of the discrete rate 
model compared to the continuous rate model. We observe that 
the discrete rate model with 8 levels performs very close to the 
continuous rate model. Moreover, as the number of the rate 
levels increases, the performance of the discrete rate model 
gets closer to that of the continuous rate model. Furthermore, 
as the number of nodes increases, the normalized schedule 
length increases, mainly due to the increasing number of 
the combinations of the links for concurrent transmissions. 
SNA-MLA performs better than SNA-MUA since it checks 
an exponential number of link combinations. 



Fig. 2. Normalized maximum total active length of SNA-MLA and 
SNA-MUA scheduling algorithms for continuous and discrete rate 
transmission models for different number of nodes at 5 nodes/m^ 
density. 

Fig.[3]shows the normalized maximum total active length of 
SNA-MLA and SNA-MUA scheduling algorithms for continu¬ 
ous and discrete rate transmission models at different network 
densities in a network of 100 nodes. Similarly, as the number 



















References 


of the rate levels increases, the performance of the discrete 
rate model gets closer to that of the continuous rate model. 
The difference between the performance of these algorithms 
decreases at low and high node densities. This is mainly due 
to the decreasing number of the combinations of the links for 
concurrent transmissions: At low node densities, the number 
of link combinations suitable for concurrent transmission is 
limited since most of the nodes are suitable for concurrent 
transmission, while causing only slight interference, being 
separated enough from each other. At high node densities, 
the number of node subsets that are separated enough from 
each other for concurrent transmission is limited due to the 
massive interference caused among each other. 



Fig. 3. Normalized maximum total active length of SNA-MLA and 
SNA-MUA scheduling algorithms for continuous and discrete rate 
transmission models at different network densities in a network of 
100 nodes. 


V. Conclusion 

In this paper, we propose a polynomial time power control 
and rate adaption algorithm for ultra-reliable M2M control 
applications employing discrete rate transmission model. The 
algorithm is based on the initialization with the vector of the 
lowest possible rates satisfying the delay constraints of the 
nodes and then checking the feasibility of the lowest possible 
number of rate allocation vectors for satisfying energy, maxi¬ 
mum power and SINR constraints. We proved the optimality 
of the proposed algorithm. The extensive simulations of 
the algorithm in combination with the previously proposed 
scheduling algorithms demonstrate that the performance of the 
discrete rate transmission model is very close to that of the 
continuous rate model with high enough discrete rate levels, 
and robust to the changes in the network size and density. In 
the future, we are planning to extend this framework for multi¬ 
hop M2M applications and next generation cellular networks. 
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